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We discuss how the the operator product expansion (OPE) can be used to derive asymptotic expressions for 
certain integrals. This yields operator matrix elements (OME's) which determine the matching conditions for 
MS parton densities across heavy flavour thresholds. Then we construct four and five-flavour densities from a 
three-flavour set via the evolution of the AP equation using LO and NLO splitting functions. 



It is well known that a s (n 2 , n/, A(n/)) in 
pQCD requires matching conditions as the scale fi 
crosses flavour matching points. At these points 
the number of light-flavours n f changes by unity 
so the QCD scale parameter A(n/) in the solution 
of the differential equation for the j3 function is 
redefined to make the running coupling contin- 
uous. When heavy quarks are included another 
scale m, the mass of the heavy quark, enters and 
the matching conditions are more complicated. 
The precise relations which need to be satisfied 
are given in [ [jj , [ || . In order lowest order pQCD 
one can choose to make the a s continuous across 
heavy flavour thresholds at /i — m. However this 
does not hold in higher order pQCD as the match- 
ing conditions in the MS scheme then contain 
non-logarithmic terms. Hence there is a disconti- 
nuity in a s at ji = m. 

Recently the analogous problem of deriving the 
two-loop matching conditions on parton densities 
as the mass factorization scale crosses the heavy 
flavour thresholds has been solved in [ ||. The 
way this was done is as follows. We examined 
the large Q 2 limit of the heavy quark coefficient 
functions which appear in NLO perturbation ex- 
pressions for heavy quark extrinsic pair produc- 
tion in deep inelastic scattering. These quanti- 
ties are functions of the virtuality of the pho- 
ton probe \J Q 2 , the mass of the heavy quark m, 
the renormalization scale fi, which is chosen equal 
to the mass factorization scale, and the partonic 
Bjorkcn scaling variable z. The number of heavy 
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D* mesons produced in deep inelastic scattering 
can be derived by convoluting these heavy (c — c) 
quark coefficient functions with appropriate com- 
binations of three-flavour light parton densities 
(u,d,s and g) and with heavy (c — c) quark frag- 
mentation functions [ 0]. Note that the heavy 
c — c pair only appears in the final state. Unfor- 
tunately we do not have analytic expressions for 
all these heavy quark coefficient functions. Some 
only exist as two-dimensional integrals over very 
complicated expressions. However there are con- 
venient tables for all of them in [ |) . 

In the limit Q 2 ^> m 2 the complicated integrals 
in the heavy quark coefficient functions reduce to 
terms with powers of ln(Q 2 //i 2 ) and \n(/j, 2 /m 2 ) 
multiplied by functions of the variable z. These 
results can be reexpressed as convolutions of 
light-mass coefficient functions C(z, Q 2 /^ 2 ) which 
contain the terms with powers in ln(Q 2 //i 2 ) and 
OME's A(z, fi 2 /m 2 ) which contain the powers in 
ln(/z 2 /m 2 ) . The way we evaluated these OME's 
is described in [ ^ so we only give an outline here. 
We wrote the heavy quark coefficient functions 
in terms of dispersion integrals for off-shell for- 
ward Compton scattering as is normally done for 
the OPE in deep inelastic scattering. We then 
changed variables to write the dispersion integral 
in terms of a variable z' which is between zero 
and unity. Next we expanded the denominator in 
a Taylor series in z'. To take the limit Q 2 3> m 2 
of the dispersion integral we add and subtract the 
same dispersion integral where we take the limit 
Q 2 3> m 2 in the integrand. This integrand con- 
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tains the OPE of the standard heavy quark (Q) 
nonsinglet and singlet operators in pQCD taken 
between states with momentum k 7 namely 

<Q(k)\O Q ^ 2 ,....^(0)\Q(k)> ■ (1) 
The heavy quark operator 

0Q, liuli2t .... lin (x) = ^(x^^D^ D^tpix) , (2) 

is a gauge invariant operator containing the heavy 
quark field tp(x) and the covariant derivative 
Dfj, = + igA^. It can be shown that the origi- 
nal integral minus the integral involving the OPE 
does not contain any mass singularities as m — » 
so it cannot depend on the heavy quark mass m 
and therefore only contains terms with powers 
of \n(Q 2 / /j, 2 ) . Hence the integrals which contain 
the evaluation of the OME's in the OPE yield 
all the terms containing powers of ln(^ 2 /m 2 ). 
This means that analytic expressions for the two- 
loop OME's with one heavy quark loop and light- 
quark or gluon incoming and outgoing states con- 
tain the information we require to extract the 
A(z,fi 2 /m 2 ). Of course the actual evaluation of 
the five OME's which exist in order a 2 requires 
the introduction of infrared and ultraviolet regu- 
lators, the use of gauge invariant operators, con- 
tractions with light-like four vectors to make the 
projections and MS renormalization. The results 
of this analysis are encapsulated in expressions 
like 

A p Q Y 2) (z,^/m 2 ) = A 1 (z)ln 2 (^/m 2 ) 
+A 2 (z)ln( f i 2 /m 2 ) + A 3 (z), (3) 

where 

A 1 (z) = C F T f [-8(l + z)hxz 

-16/(3z)-4 + 4z + 16z 2 /3], (4) 

A 2 (z) = C F T f [-8(1 + z)]v?z 

+ (8 + 40z + 64z 2 /3) In 2 

+ 160/(9z)-16 + 48z-448z 2 /9], (5) 

A 3 (z) = C F T f {(l + z)[32S h2 (l-z) 
+ 161nzLi 2 (I - z) 
-16C(2)lnz-4/31n 3 z] 



+(32/(32) + 8 - 8z - 32z 2 /3)Li 2 (l - z) 

+(-32/(32) - 8 + 8z + 32z 2 /3)C(2) 

+ (2 + 10z + 16z 2 /3) ln 2 z 

-(56/3 + 88z/3 + 448z 2 /9) In 2 - 448/(27z) 

-4/3- I24z/3+ 1600z 2 /27}. (6) 

The tilde indicates that an overall factor of n/ 
has been extracted from the function and the (2) 
in the superscript means this is the second or- 
der term in an expansion in a s — a s /(47r). The 
five functions A^ q (z, /i 2 /to 2 ), where PS de- 
notes pure singlet under the flavour group (i.e., no 
non-singlet projection exists), Aq^\z, // 2 /m 2 ), 
where S denotes singlet under the flavour 
group, A S gg {2 ^(z,fi 2 /m 2 ), A S gq (2 ^(z, fi 2 /m 2 ), and 
A^ g Q 2 \z, ^ 2 /m 2 ) where NS denotes non-singlet 
under the flavour group, which exist in order a 2 
pQCD are given in [ |). Alternative discussions 
of their derivation and use are given in [§. Note 
that they contain nonlogarithmic terms such as 
A3 (2) in Eq.(6) in order a 2 so there is no scale 
fi where we can make them all vanish. Since we 
know the four-flavour light mass coefficient func- 
tions C(z, Q 2 j ' [i 2 ] in order a 2 [| we can analyt- 
ically evaluate the convolutions with the appro- 
priate A's to obtain asymptotic expressions for 
the heavy quark coefficient functions. They were 
given in [ |9|. As far as this workshop is con- 
cerned we would like to point out that this "in- 
verse mass factorization method" is an elegant 
use of the OPE to obtain asymptotic expansions 
of integrals. 

Normally parton densities are fitted to specific 
functions of a; at a scale /x and the AP equations 
then govern the evolution of these densities to 
other scales. Suppose one begins with a three- 
flavour set containing densities for u,d,s quarks 
and the gluon g. Then the above results allow one 
to define four-flavour parton densities at scales 
fi > m c from the input set of three-flavour den- 
sities in fixed-order perturbation theory (FOPT). 
Let the <g> symbol denote the convolution integral 
f®9 = ! f(x/y)g(y)dy/y, where x < y < 1, then 
we define the charm density 

fc+c{n f + 1, fJ, 2 ) = 
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Figure 1. £Cnnlo(5, x, [i 2 ) in the range 10 5 < 
x < 1 for fi 2 = 20.25, 25, 30, 40 and 100 in units 
of (GcV/c 2 ) 2 . 
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a s {n f ^ 2 )A% g (^)®f g {n f ^ 2 ) 
WAn f ^ 2 )[Al\(^)®f q {n h f) 



+4 ff (£ ! )®/K% ) M 2 ) 

the singlet gluon density 

/ g s (n / + l, M 2 )=/ s s (n/5/i2) 

u 2 

+a s (n f ^ 2 )A s ggQ {^)®f g {n f ^ 2 ) 
m c 

r a* 2 

+a 2 (n/, /i 2 ) [ A gq ,Q(^2 ) ® fl( n f^ 2 ) > 



(7) 



(8) 



and the light mass quark densities 
4+fc( n / + !'M 2 ) = fk+k(n f ,^ 2 ) 

2 

+a 2 (n / , A1 2 )^ s Q (^) ®/ fe+ fe(n /;M 2 ), (9) 

V 771 c / 

for = 3 and m 2 < /i 2 < m 2 . Note that we 
have suppressed the a; dependence to make the 
notation more compact. These expressions were 
used in [ [?J to construct a variable flavour number 



Figure 2. Same as Fig.l for the NLO results from 
MRST98 set 1. 

scheme (VFNS) for the heavy quark contributions 
to the deep inelastic structure functions. 

Note however that the above procedure 
does not resum the potentially large terms in 
ln(^ 2 /m 2 ) which are explicitly left in the par- 
ton densities. To do this we need to evolve the 
above densities via the AP equation rather than 
using FOPT. This is new work in [ [w| using three- 
flavour densities at small scales from [ |ll| . The 
latter LO and NLO densities are started at very 
small scales fio below the mass of the charmed 
quark. Hence three flavor evolution proceeds 
from the initial /Xq to the scale [i 2 = m 2 = 1.96 
(GeV/c ) 2 . In this region a s is large so we had to 
be very careful to get numerically accurate solu- 
tions of the evolution equation. Fortunately there 
are standard inputs and tables in [ [l2| with which 
we could compare the parton densities from our 
evolution code. We chose the matching scale \x 
at the mass of the charm quark m c so that all 
the ln(/i 2 /m 2 ) terms in the OME's vanish at this 
point leaving only the nonlogarithmic pieces in 
the order a 2 OME's to contribute to the right- 
hand-sides of Eqs. (7), (8) and (9). Note that 
the LO and NLO charm densities vanish at the 
scale fi — m c since 

4f(z, ^/rn 2 ) = 4T/(z 2 + (1 - z f) ln^ 2 /m 2 ) , 
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Figure 3. Same as Fig.l for the NLO results from 
CTEQ5HQ. 

does not have a non-logarithmic term. The 
NNLO charm density starts off with a finite in- 
dependent shape in order a 2 determined by 

fc+c(n f + l,m%) = 



a 2 s (nf,m 2 c ) 



+A S Qg (l) 



A™{l)®f q {n f ,ml) 



(11) 



with rif = 3. Hence the OME's provide the 
boundary condition for the evolution of the 
(massless) charm density. Also note that we or- 
dered the terms on the right-hand-side of Eq.(ll) 
in powers of a s so that the result contains a prod- 
uct of NLO OME's and LO parton densities, al- 
though this is not evident here. The result is 
then strictly order a 2 and should be multiplied 
by order a° coefficient functions when forming the 
zero-mass variable flavour number scheme (ZM- 
VFNS) charm density contribution to the deep 
inelastic structure functions. 

The four-flavour gluon density is also generated 
at the matching point in the same way. At /i = 
m c we define 

f^n f + l,m 2 c ) = f^(n f ,m 2 c ) 



+a 2 s (n f ,m 2 c ) 



A S gq , Q (l)®fq(n f ,m 2 c ) 



Figure 4. x6nnlo(5, x, [i 2 ) in the range 10 5 < 
x < 1 for [i 2 = 20.25, 25, 30, 40 and 100 in units 
of (GcV/c 2 ) 2 . 

The four-flavor light quark (u,d,s) densities are 
generated using 

fk+k( n f + l' m c) = /fc+fc( n /> m c) 

+a 2 s (n f , m 2 c )A™ Q ( 1 ) ® f k+ - k (n f , m 2 c ) . (13) 

The total four-flavor singlet quark density follows 
from the sum of Eqs.(ll) and (13). 

Next the resulting four-flavor densities are 
evolved from their boundary values using the 
four-flavor evolution kernels in the AP equations 
in either LO or NLO up to the scale /j 2 — 20.25 
(GeV/c 2 ) 2 . The bottom quark density is then 
generated at this point using 

fb+b(nf + l,ml) = 

a 2 (n f , m 2 ) [i^ (1) ® /f (n f , m 2 ) 

+ig(l)®/ 9 s (n / ,m 2 )" , 



(14) 



+4 g , Q {l)®f g {n f ,ml) 



(12) 



and the five-flavour gluon and light quark den- 
sities (which now include charm) are generated 
using Eqs. (12) and (13) with n/ = 4 and re- 
placing m 2 by m\. Therefore only the nonloga- 
rithmic terms in the order a 2 OME's contribute 
to the matching conditions on the bottom quark 
density. Then all the densities are evolved up to 
higher fj 2 as a five-flavor set with either LO or 
NLO splitting functions. 
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Fig. 5 
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Figure 5. Same as Fig. 4 for the NLO results from 
MRST98 set 1. 

The above formulae and their evolution with 
LO and NLO splitting functions have been imple- 
mented in a C++ computer code [ [to] to yield the 
CS parton density set. They were used in the con- 
struction of two VFNS for the charm quark con- 
tribution to the deep inelastic structure functions 
in [ . Note that approximate expressions for 
the three loop splitting functions are now avail- 
able in [ Q . When NNLO parton densities are 
available from fits to experimental data we can 
incorporate them into our computer program. 

As an illustration we would like to compare the 
charm and bottom quark densities in the CS [ |To) , 
MRST98 [ H| and CTEQ5 [ @ sets. The lat- 
ter two sets work with order a s matching condi- 
tions so the parton densities are continuous across 
heavy flavour thresholds. The MRST98 sets use 
a procedure proposed in [ 0, while the CTEQ5 
sets use the different ACOT procedure in [ |l8[ . 
Here we show the five-flavor densities. In the CS 
set they start at fi 2 =m\ = 20.25 GeV 2 . At this 
scale the charm densities in the CS, MRST98 (set 
1) and CTEQ5HQ sets are shown in Figs.1,2,3 
respectively. Since the CS charm density starts 
off negative for small x at fi 2 = m 2 = 1.96 GeV 2 
(see the plots in [ ) it is smaller than the corre- 
sponding CTEQ5HQ density. At larger /i 2 all the 
CS curves in Fig.l are below those for CTEQ5HQ 
in Fig. 3 although the differences are small. In 
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Figure 6. Same as Fig. 4 for the NLO results from 
CTEQ5HQ. 



general the CS c-quark densities are more equal 
to those in the MRST98 (set 1) in Fig. 2. 

At the matching point fi 2 = 20.25 GeV 2 the fa- 
quark density also starts off negative at small x as 
can be seen in Fig. 4, which is a consequence of the 
explicit form of the OME's in [ |. At 0(a 2 ) the 
nonlogarithmic terms do not vanish at the match- 
ing point and yield a finite function in x, which 
is the boundary value for the evolution of the b- 
quark density. This negative start slows down 
the evolution of the b-quark density at small x as 
the scale n 2 increases. Hence the CS densities at 
small x in Fig. 4 are smaller than the MRST98 (set 
1) densities in Fig. 5 and the CTEQ5HQ densities 
in Fig. 6 at the same values of /j, 2 . The differences 
between the sets are still small, of the order of 
five percent at small x and large /i 2 . This will 
lead to differences in cross sections for processes 
involving incoming b-quarks at the Tevatron. 

We suspect that the differences between these 
results for the c and b-quark densities are pri- 
marily due to the different gluon densities in the 
three sets rather to than the effects of the differ- 
ent boundary conditions. This could be checked 
theoretically if both LO and NLO three-flavor 
sets were provided by MRST and CTEQ at small 
scales. We note that CS uses the GRV98 LO 
and NLO gluon densities, which are rather steep 
in x and generally larger than the latter sets 
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at the same values of /j, 2 . Since the discontinu- 
ous boundary conditions suppress the charm and 
bottom densities at small x, they enhance the 
gluon densities in this same region (in order that 
the momentum sum rules are satisfied). Hence 
the GRV98 three flavour gluon densities and the 
CS four and five flavor gluon densities are gen- 
erally larger than those in MRST98 (set 1) and 
CTEQ5HQ. Unfortunately experimental data are 
not yet precise enough to decide which set is the 
best one. 
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